Directed d-mer diffusion describing Kardar-Parisi-Zhang type of surface growth 
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We show that d + 1-dimensional surface growth models can be mapped onto driven lattice gases of 
d-mers. The continuous surface growth corresponds to one dimensional drift of d-mers perpendicular 
to the (d — l)-dimensional "plane" spanned by the d-mers. This facilitates efficient, bit-coded algo- 
rithms with generalized Kawasaki dynamics of spins. Our simulations in d — 2, 3, 4, 5 dimensions 
provide scaling exponent estimates on much larger system sizes and simulations times published so 
far, where the effective growth exponent exhibits an increase. We provide evidence for the agree- 
ment with field theoretical predictions of the Kardar-Parisi-Zhang universality class and numerical 
results. We show that the (2 + 1) -dimensional exponents conciliate with the values suggested by 
Lassig within error margin, for the largest system sizes studied here, but we can't support his 
predictions for (3 + 1) d numerically. 
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One of the simplest nonlinear stochastic differential 
equation set up by Kardar, Parisi and Zhang (KPZ) jlj 
describes the dynamics of growth processes in the ther- 
modynamic limit. It specifies the evolution of the height 
function /i(x, t) in the d dimensional space 

d t /i(x, t)=v + o-V 2 /i(x, t) + A(V/i(x, t)f + ?7(x, t) . (1) 

Here v and A are the amplitudes of the mean and local 
growth velocity, a is a smoothing surface tension coeffi- 
cient and r\ roughens the surface by a zero-average Gaus- 
sian noise field exhibiting the variance (^(x, t)r)(x.' , t')) = 
2D5 d (x — x')(t — t'). The notation D is used for the noise 
amplitude and () means the distribution average. 

The KPZ equation was inspired in part by the the 
stochastic Burgers equation [2J, which belongs to the 
same universality class @, and it became the subject of 
many theoretical studies Besides, it models other 

important physical phenomena such as directed polymers 
Q , randomly stirred fluid @ , dissipative transport [1, Q , 
and the magnetic flux lines in superconductors [Tcl | . The 
equation is solvable in (1 + 1) d [11|, but in higher dimen- 
sions approximations are available only. As the result 
of the competition of roughening and smoothing terms, 
models described by the KPZ equation exhibit a rough- 
ening phase transition between a weak-coupling regime 
(A < A c ), governed by the Edwards- Wilkinson (EW) 
fixed point at A = pjj], and a strong coupling phase. 
The strong coupling fixed point is inaccessible by pertur- 
bative renormalization group (RG) method. Therefore, 
the KPZ phase space has been the subject of controver- 
sies and the value of the upper critical dimension has 
been debated for a long time. 

Using a directed polymer representation, the validity of 
a scaling hypothesis [13[ and the two-loop RG calculation 
for d > 2 [14j was confirmed and extended to all orders 
in d = (2 + e) These results provided an argument 
for an upper critical dimension d c — 4 of the roughening 
transition, but the strong-coupling rough phase is not 
accessible by perturbation theory. Above d = 1 the scal- 



ing behavior in the rough phase has been very controver- 
sial, diverse values for the scaling exponents were claimed 
(16h18| . In particular, assuming that height correlations 
exhibit no multiscaling and satisfy an operator product 
expansion, exact field-theoretic methods lead to ratio- 
nal number growth values in two and three dimensions 
[l~7j . Some theoretical approaches predict that d c = 4 is 
an upper critical dimension of the rough phase [l9l - [22l |. 
Recently, a non-perturbative RG study has been able to 
describe the strong coupling fixed point and has provided 
indications for a possible qualitative change of the criti- 
cal behavior around d = 4 123 1 . This is in contradiction 
with the numerical results [24| - l27| , which predict the lack 
of an upper critical dimension. 

Mapping of surface growth onto reaction-diffusion sys- 
tem allow effective numerical simulations [HI, [2t| . As a 
gen eralization of the 1 + 1 dimensional roof-top model 
3C . Hl| and the 2 + 1 dimensional octahedron model 
32j, we consider the deposition and removal processes 



of higher dimensional objects on d > 2 dimensional sur- 
faces. We remind that in 1 + 1 dimensions a continuous 
surface line having no overhangs can be approximated 
by 45 degree up/down slope elements (after appropriate 
length rescaling). A process with KPZ scaling can be 
realized by deposition at local minima or removal of lo- 
cal maxima (roof-top) . If we associate the up slopes with 
'particles' and down slopes with 'holes' of the base lattice 
(see Fig. [TJa), the adsorption/desorption corresponds to 
the asymmetric exclusion process (ASEP) 33]. This is a 
lattice gas [H| [35|, where particles can hop on adjacent 
sites with asymmetric rates and hard-core exclusion. Its 
behavior is well known, and variations of ASEP (disor- 
der, interactions ... etc.) correspond to variations of 1 + 1 
dimensional KPZ growth models. 

We have extended the roof-top construction to (2 + 1) 
dimensions 32] by the introduction of octahedra having 
four slopes. The up edges in the x or y directions can be 
represented by '+1', while the down ones by '—1', and 
a surface element update is a generalized (Kawasaki) ex- 
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change (Eq. (3) of [32]). The translation of up edges to 
'particles' and the down ones to 'holes' of the base lattice 
maps particle deposition/removal processes onto two si- 
multaneous particle moves (one in the x and one in the 
y direction) . One can also consider it as a dimer move in 
the bisectrix direction of x and y (see Fig. [TJb) . There- 
fore, the (2 + 1) surface dynamics can be mapped onto 
a "two-dimensional ASEP" of oriented dimers exhibiting 
hard-core exclusion. The asymmetric drift corresponds 
to an evolving surface exhibiting KPZ scaling, while the 
symmetric dimer diffusion is related to the EW behavior. 

Now we proceed with this kind of construction, con- 
sidering the discrete slope variables in higher dimensions 
and generalize the simultaneous +1 <H> — 1 (Kawasaki) 
exchange rule of them (Eq.(3) of [HJ ) to d-dimensional 
updates 



(2) 



with probability p for attachment and probability q for 
detachment (see Fig.[TJc for the 3d case). It is well known 
[i| that the surface evolution of the deterministic KPZ 
growth are described also by the Burgers equation [2] for 
growth velocities v(x,t) in the surface normal obeying 

<9 t v(x, t) = ctV 2 v(x, t) + Av(x, t) Vv(x, t) (3) 

due to the transformation v(x, t) — Vh(x,t). 

In the forthcoming part we will prove that our micro- 
scopic model for <i-mers in the continuum limit can be 
mapped onto the anisotropic version of Eq. ([3]), similarly 
as shown in lower dimensions [3ll l32j . The derivation is 
based on the formulation of the reduction of possible up- 
dates. Our surface model is represented by the discrete 
derivative elements: 5 X , S y , S z ... (e ±1) at every lattice 
points. A generalized Kawasaki update @ is defined by 
a matrix 



/ S x (i - fc,...) 
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Uy{l,j - l,fc, ...) 

6 z (i,j,k-l,...) S z (i,j,k,...) 
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In d dimensions we define vectors of the slopes, the 
columns of (U]), analogously to one and two dimensions: 
0i,j,k,~ = ~ 1, J, fc, —),8 y (i,j - l,fc, ...), ...) , around 
the lattice point, which we select for deposition/removal 
update and set up a microscopic master equation 

d t P({W},t) = <j,k,...(M)P({v},t) 

i,j,k,... 

E w i,i,k,.XM)P({°ht) (5) 

i,j,k,... 

with the probability distribution P({a},t). Here the 
prime index denotes the state of a following the update 



([2]). The transition probability of as can be expressed 
as 

w i,j,k,...({?}) = — O r i+lJ4-l,fc+l,...O r iJ,fc,... (6) 

+ ^(&i+l,j + l,k+l,... ~Vi,j,k,...)] j 

with A = 2p/(p + q) — 1 parametrization, which formally 
looks like the Kawasaki exchange probability in Id, ex- 
cept the factor A, which is necessary to avoid surface 
discontinuity creation in higher dimensions. This means 
that we update the slope configurations only if the values 
of all coordinates of the vector W are identical as shown 
by @. One can allow formally these updates via the 
expression 

A = l/2 d+1 det[(a t ^ K .+Ca h3 ^..) (7) 

x (&i+l,j+l,k+l,.. + C^<+l,3+l,A+l,..)iJ i 

where / and C are the unity and the cyclic permuta- 
tion matrices respectively. The matrix C shifts each 
coordinate value to the next index value. Thus for W- 
s with mixed coordinate values, the vectors (cFij fc + 

Ca%,3,k,..) = fc or (o^+ij+i.fc+i,.. + (?o r t+i,i+i,fc+i 1 ..) = 
possess zero elements. Therefore the determinant of kk'I, 
being the product of the diagonal elements, is zero in case 
of mixed coordinates and A = l/2 d+1 in case of equal co- 
ordinates. 

For example a d = 3 update is prohibited when the 
slope vector is ~o = (1,1,-1), because k has one coordi- 
nate value of zero 
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By calling '+l'-s as particles and the '— l'-s as holes 
of the base lattice, their synchronous update can be con- 
sidered to be a single step motion of an oriented cZ-mer 
in the bisectrix direction of the x, y, z, ... coordinate 
axes. Thus <i-mers follow one-dimensional kinetics, de- 
scribed by Kawasaki exchanges To obtain a one- 
to-one mapping we update neighborhoods of the lattice 
points denoted by the green dots of Fig. [T] 

To derive Eq. ([3]) first we have to average over the slope 
vectors 



5>P({7f},t). 

{5} 



(8) 



By calculating its time derivative using the master equa- 
tion ([5]) and the transition probabilities ([6]) 



{a} i,3,k,... 

a E «^-, fc ,...({<7})P({<7},t)] , 

i,j,k,... 



(9) 



in which we filter out vectors of non-equal coordinates 
([7]) (thus Wi.j^k,... is nonzero only if <cFi,j,... i= o'i+i 
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FIG. 1: (Color online) Mapping of one (a), two (b), three 
(c) dimensional surface growth (discrete up/down derivatives) 
onto 1, 2, 3-dimensional oriented motion of d-mers. In general 
a deposition event corresponds to a d-Kawasaki exchange, or 
a d-mer hopping in the bisectrix direction of axes. Up/down 
slopes between lattice points (green dots) are denoted by 
full/empty circles, corresponding to particles/holes on the 
projection space. The surface continuity translates into par- 
ticle exclusion, which makes the model nontrivial in higher 
dimensions. 



we can obtain 

2dt(vi,j,k,...) = (o^-ij-i,...) - (10) 

+ (^i+lj+l,...) + A(<7j )J -(<7j-|_i ) j-_|_i ) ... — , 



analogously to one dimension [31] . Here one can see the 
discrete first and second differentials of Wi.j,k,... corre- 
sponding to the operators of Eq. ^ in the bisectrix direc- 
tion of the axes. These differentials are one-dimensional, 
because the cZ-mer dynamics is one-dimensional. In prin- 
ciple one could derive a set of coupled Burgers equations 
for the particles in each direction in an isotropic way in 
accordance with isotropic surface model, but the coordi- 
nated movements reduce the dimensionality and we can 
map onto an anisotropic equation of d-mers. 

Making a continuum limit in each direction and taking 
into account the relation of height and slope variables 
(v(x, t) — V/i(x, t)), we can arrive to the deterministic 
KPZ equation. The nonlinear term vanishes for p = q 
(A = 0). The sign of the coefficient A of the nonlinear 
term can be interpreted as follows: For p > q positive 
non-linearity (positive excess velocity) it is a consequence 
of growth with voids. 

Since this derivation was applied just for the first one 
in the hierarchy of equations for correlation functions it 
does not prove the equivalence to the stochastic KPZ. 
Furthermore, the form of the noise term, which was not 
considered in our derivation, may also introduce differ- 
ences. Although our surface model is spatially isotropic, 
we can map it onto a one-dimensional Burgers equation 
(of extended objects), therefore anisotropic scaling be- 
havior might be expected. However, by going into the 
continuum description the hard-core exclusions necessary 
to provide continuous surfaces is lost, and the resulting 
equation looks trivial. 



Here we investigate by numerical simulations this iso- 
topic surface growth model via the one-dimensional di- 
rected migration of eZ-mers in the d-dimensional space. 
We have developed bit-coded algorithms for the updates 
([2]) and run it with p = 1, such that randomness comes 
from the site selection only. Therefore it is important to 
use a very good random number generator, which pro- 
vides uniformly distributed numbers with high resolu- 
tion. Otherwise we would realize a process with quenched 
disorder, which for KPZ belongs to a different universal- 
ity class (s ee 12910 . We used the latest Mersenne- Twister 
generator [36| in general, which has very good statisti- 
cal properties and which is very fast, especially in the 
SSE2 registers. We tested our results using other ran- 
dom number generators as well. In practice each update 
site can be characterized by the 2 different local slope 
configurations. However, due to the surface continuity 
we need only a few bits of a world (1 byte for d = 2, 3, 4 
and two bytes for d = 5) for this purpose. This allows 
an efficient storage management in the computer mem- 
ory and permits simulations of larger system sizes. The 
updates can be performed by logical operations, either 
on multiple samples at once, or on multiple (not overlap- 
ping) sites at once. Our bit-coded algorithm proved to be 
~ 40 times faster than the conventional FORTRAN 90 
code we started with. It is important to note that this 
stochastic cellular automaton like representation of the 
surface growth opens the possibility for an implementa- 
tion on extremely fast graphic cards with massively par- 
allel processors. Furthermore the construction permits 
the extension of the mapping for more complex surface 
reactions (37j . 

We performed dynamical simulations by starting from 
stripe ordered particle distributions. This corresponds 
to a flat surface with a small intrinsic width. The 
considered lattices gases had the maximum linear sizes 



L max = 2 15 ,2 10 ,2 8 ,2 6 for d = 2,3,4,5 dimensions, re- 
spectively and periodic boundary conditions were ap- 
plied. A single step of the lattice gas algorithm com- 
prises a random site selection and in case of an appropri- 
ate neighborhood configuration a p = 1 Kawasaki <i-mer 
update ([2]). The time is incremented by 1/L d in units of 
Monte Carlo steps (MCs). Throughout the paper we will 
use this unit of time. 

We could exceed by magnitudes of order all previous 
numerical system sizes and simulation times. For ex- 
ample the largest five-dimensional simulations were done 
for L = 30 and t m = 230 MCs [25]. Our L = 64 simu- 
lations, where we have the good bulk/surface ratio: 5.4, 
required 2 GB memory size and a couple of weeks for a 
single realization up to t max = 5000 MCs. Similarly, the 
largest sized simulations in d = 2 for L = 1 1520 sys- 
tem could achieve t max = I0 4 MCs (38|. Our largest 
L = 32768 sized simulations reached t max = 44600 MCs. 
The longest runs for L = 4096 passed the saturation at 
t ~ 4 x I0 5 MCs and the samples were followed up to 
t rnax = I0 6 MCs. 

We run the these lattice gas simulations for 10 — 1 000 
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independent realizations for each dimension and size con- 
sidered, and calculated h x y (t) and the second moment 



W(L,t) 
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(11) 

from the height differences at certain sampling times. 
The growth is expected to follow the Family- Vicsek scal- 
ing |39j asymptotically, but due to the corrections it can 
be described by a power series 



W(t,L -»• oo) = bt f3 (l + b t^° + hit* 1 ...) 



(12) 



with the surface growth exponent /?. For finite system, 
when the correlation length exceeds L, the growth crosses 
over to a saturation with the scaling law 



W(t -+oc,L) = aL a (l + a L w ° + ai L Ul ...) 



(13) 



characterized by the roughness exponent a. In our case 
the intrinsic width of the initial state, which is repre- 
sented by a zig-zag surface of width 1/2 (see Fig FT]), 
results in a constant correction term. Thus we have 
&o = 1/2, 0o = —j3 and ciq = 1/2, luq = —a. Dur- 
ing our scaling analysis we dropped this contribution by 
subtracting W 2 (0) = 1/4 from the raw data and consider 
the next leading order correction as leading one. Further- 
more we disregarded the initial time region t < t ~ 50, 
when basically an uncorrelated random deposition oc- 
curs. The dynamical exponent z can be expressed by the 
ratio z — a/ /3 and in case of the Galilean invariance of 
an isotropic KPZ equation the z = 2 — a relation should 
also hold. 

Besides the extensive simulations we have performed 
careful correction to scaling analysis by calculating the 
local slopes of the exponents. The effective exponent of 
the surface growth can be estimated similarly as in case 
of other scaling laws [2{|, as the discretized, logarithmic 
derivative of (fTTI) 



a eff (L) 



In W(t ->ao,L)- In W(t -)• oo, L') 
ln{L) - ln(L') 



(14) 



It was determined numerically for different discretiza- 
tions: t/t' — 2,3, and we tried to fit it with the leading- 
order correction ansatz, which can easily be deduced from 
p| (see 13 or [H) 



/3 e// (i)= y 8 + 6 1 ^ 



(15) 



for t > to and before the saturation region. In other 
cases, such as ballistic deposition, which has a large 
unknown intrinsic width one can use another effective 
roughness exponent definition introduced in [40(. 

We tested our method with the one- dimensional, ex- 
actly known case. Simulations were run on L = 5 x 10 D 
sized system up to t max = 16666 MCs for 40 indepen- 
dent realizations. We determined the effective exponents 
j3 e ff(t), which approaches /3 — 1/3 from below, in a 
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FIG. 2: (Color online) Effective exponents of the growth of 
the d = 1 dimensional model. The solid line corresponds to 
the simulation result. The dashed lines shows a fitting with 
the form (Tl5)) . 



perfect agreement with the leading-order correction form 
(see Figure^]). The fitting with (TT5|) on the local slopes 
data resulted in fj = 0.333(5) and (j>i = —0.53. 

Similarly to the time dependence we can analyze the 
size dependence following the saturation by determining 
the effective exponent of the roughness, which can be 
defined as the logarithmic derivative of (fTTj) 



a eff(L) 



In W(t oo, L) - In W{t -> oo, V) 
ln(L) - ln(L') 



(16) 



The finite size scaling was done for systems of linear sizes 
in between L max (discussed earlier) and L min , which was 
2 6 for 2d, 2 5 for 3d, 2 4 for 4d and 2 3 for 5d, respectively. 
To handle the boundary conditions effectively, system 
sizes of power of 2 were simulated. To get the asymptotic 
values we took into account all effective exponent points 
shown on Fig. U and applied a leading order, linear fit- 
ting. The error margins of exponents are estimated from 
the error-bars of Fig. 01 This method gives a better es- 
timate for the asymptotic values than just a least-square 
fitting on the data points, which completely disregards 
corrections to scaling. We also calculated rough, but in- 
dependent estimates for z by measuring the relaxation 
time, i.e. the time needed to reach 90% of the saturation 
value. The asymptotic value is extrapolated by a linear 
fitting: z e ff(t) = z + C\/L. 

In two dimensions we estimated the growth exponent 
in the largest system sizes considered (L — 2 15 , 2 14 , 2 13 ) 
(see Fig. 0). Fitting in the 50 < t < 44600 time 
window with the form (fTB")) resulted in b\ = 0.83 and 
/3 = 0.245(5), which is somewhat bigger than what was 
obtained by the largest known (L = 11520) sized sim- 
ulations: f3 = 0.240(1) [11], and all other previous nu- 
merical estimates including ours [H, EH SH ■ This value 
conciliates with the f3 — 1/4 RG exponent of 1_7[. One 
can obtain this value by the late time behavior of effec- 
tive exponent, which has not been seen before, because 
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FIG. 3: (Color online) Effective exponents of the growth of 
the d = 2,3,4,5 dimensional model (top to bottom). Solid 
lines correspond to the largest system L — 2 15 , L — 2 10 , 
L = 2 8 and L — 2 6 respectively. Dashed lines show our 
results of smaller sizes: L = 2 10 , L = 2 s , L = 2 5 and L = 2 5 
respectively, where saturation sets in earlier causing a cutoff 
in the scaling. 

finite size effects have screened it. On the graph one 
can see strong oscillations for L = 2 10 and intermediate 
times, which are damped before saturation. In the one- 
dimensional ASEP model such oscillations are shown to 
be the consequence of density fluctuations being trans- 
ported through a finite system by kinematic waves (42j . 
One can speculate that the slight final increase of j3 e f / for 
the largest system sizes is just a fluctuation or oscillation 
effect, but we could not eliminate this overall tendency 
by increasing the statistics. Although the statistical fluc- 
tuations grow dramatically, as t — > oo the increase of the 
mean value is observable for each size L > 2 . Our error- 
bar of P reflects this uncertainty. The width saturation 
values have been investigated for L = 2 6 , 2 7 , 2 12 . We 
took into account the leading order correction to-scaling 
by the following Ansatz 

a e ff(L) = a + aicuiL^ 1 , (17) 

but due to the larger error-bars we restricted it to a linear 
approximation: ui = — 1. The local slopes of the steady 
state values a e //(l/L) and of z e ff(l/L) are shown on 
Fig.H This provides a = 0.395(5) and oi = 2.02 for the 
roughness and z e ff = 1.58(10), with the linear coefficient 
Ci = 1.83 for the dynamical exponent. This roughness 
exponent is in agreement with RG value (l7j . and some- 
what bigger than the existing figures a = 0.393(3) [26[ 
for L < 1024 and a = 0.385(5) 38] for L < 128. 

In three dimensions the local slope analysis for L = 2 10 
results in bi — 0.1 and j3 — 0.184(5) agreeing with the nu- 
merical results from the literature: /3 — 0.180(2) [251138$. 
P = 0.186(1) 1261 . But our estimate is much higher than 
P = 0.168(3) [4l[ (based on L < 200 sized simulations) 
and P = 1/6 predicted by RG [ItJ • For the saturation we 
obtained a\ — 1.40 and a = 0.29(1), matching a — 0.29 



FIG. 4: (Color online) Local slopes of the finite size scaling of 
the saturation width in d = 2, 3, 4, 5 dimensions (top to bot- 
tom). Inset: Effective exponents of the characteristic times. 



of 411 and in marginal agreement with a — 0.3135(15) of 
[26l | and a — 0.308(2) of [25|. The direct z e ff measure- 
ment exhibits a strong correction to scaling: z = 1.60(1) 
(ci = 1.10) and one cannot differentiate it from the 2 + 1 
dimensional results within the error margins. 

In four spatial dimensions our best fit for the growth 
exponent is b\ = 1.08 and (5 — 0.15(1). In the liter- 
ature P = 0.16(1) 13 and p = 0.146(1) [H values 
are reported. For the width saturation values the lin- 
ear fitting results in a = 0.245(5) with a\ = 0.07. This 
compares with the literature values a = 0.245(1) [25[ 
and a = 0.255(5) [26|. The z e tf seems to converge to 
z = 1.91(10) (with ci = —0.64) but the fluctuations are 
very strong and we could not reach saturation for sizes 
larger than L = 128. Going further by a factor of two in 
system sizes would require simulations with 8GB memory 
and very long CPU times. Our results do not support the 
field theoretical prediction of d c = 4, because we don't 
observe the disappearance of power-law growth. 

In five dimensions the local slopes suggest b\ = 0.134 
and p = 0.115(5) in agreement with p = 0.11(1) [25| re- 
ported for smaller sizes. One can find strong oscillations 
before the saturation regime. Again these are due to 
kinematic transport waves in finite system. Initially for 
L = 64 we saw a definite increase in /3 e // as 1/t < 0.005 
before the saturation, but this proved to be an artifact of 
the MT random number generator. When we used dif- 
ferent, pseudo-random number generators: drand48 [43[ 
or random() of language C, the growth tendency for very 
late times was much weaker. We think that the site se- 
lection, the only source of randomness in case of p = 1 
might not be completely uniform among the 2 30 possible 
places. To confirm this we repeated the 5d simulations 
using p — 0.9 with the MT generator and found agree- 
ment with the results using drand48. For the saturation 
we estimate a = 0.22(1) with a\ — 0.08 and z — 1.95(15) 
with ci = —0.55. 



6 



In conclusion we have shown that the mapping of 
a KPZ surface growth model onto driven lattice gases 
(DLG) can be extended to higher dimensions. Although 
the growth of the surfaces exhibits the spatial symme- 
try of the underlying lattice, one can map it onto an 
anisotropic DLG of more complex objects. The coarse 
grained, continuum description of these <i-mers is an 
anisotropic Burgers equation. Still the DLG model is 
non-trivial, because it is just an oriented drift of d-mers 
with hard-core exclusions. The topological constraint 
is the consequence of the required surface continuity by 
the mapping. In two dimensions we confirmed (37j that 
the probability distribution P(W 2 ) matches the univer- 
sal scaling function determined for another KPZ model 
[27j . We presented effective bit-coding simulations and 
high precision results for the exponents a, ft and z inde- 
pendently (see Table I.). The sensitive local slope anal- 
ysis provides numerical agreement with former simula- 
tion results, but for larger sizes, which have not been 
investigated so far, we see aslight growing tendency in 
the f) e ff exponents in all dimensions. For d = 2 our 
results marginally overlap with the /3 = 1/4 value sug- 
gested some time ago by RG. The change towards a triv- 
ial behavior in higher dimensions in the DLG language 



would mean the disappearance of the topological con- 
straints among the extended cJ-mer objects as they could 
follow a simple ASEP dynamics of point particles. This 
will be the target of further studies using massively par- 
allel algorithms on graphic cards. We hope that we will 
be able to obtain a firm estimate for the upper critical 
dimension using extrapolation techniques. 

TABLE I: Independent growth exponent estimates of the ri- 
mer model in different dimensions 



(/ 


a 


P 


z 


2 


0.395(5) 


0.245(5) 


1.58(10) 


3 


0.29(1) 


0.184(5) 


1.60(10) 


4 


0.245(5) 


0.15(1) 


1.91(10) 


5 


0.22(1) 


0.115(5) 


1.95(15) 
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